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Abstract 


This paper reports on the development of a compact radiation monitor/dosimeter, the LIULIN-3M, 
and on extended measurements conducted on the ground and on commercial aircraft on domestic 
and international flights. 
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I. INTRODUCTION 

The LIULIN-3M is a “Low-LET Radiation 
Spectrometer” (LoLRS). It evolved from an 
international cooperative project by a group 
of Bulgarian, Russian, German, and Ameri- 
can Scientists. The radiometer is a low- 
power, small size, lightweight, and low-cost 
instrument composed of a solid state detec- 
tor (SSD) with supporting electronics that 
enable it to operate as a pulse height ana- 
lyzer of energy deposited in the detector, 
and to obtain from these measurements the 
total dose or the dose rate produced by 
charged particles. A flash memory allows 
self-storage of data during flights and post 
flight retrieval [1 ,2,3,4], 

n. DESCRIPTION 



H— 


DC/DC converter 
SSD 62 V Bias 


A-225 Ampl 
240 mV/MeV 

X 


K- 


External Test 
Generator 


Discriminator Threshold: 38.7 mV 

— x 


Peak Hold Detector 

— x 


ADC; 256 Steps @13 mV 

I 


Power Supply 
2 x 3.6 V, 13 Ah 
(LSH20 type) 


SRAM 512 BYTES U-J 


Microcontroller 

PIC16C74 


Parallel 

Interface 


JL 


k-4 Microcontroller PIC16C74 


Flash Memory 
1 MB 


Figure 1. LIULIN Block Diagram 


The basic detector is a Silicon- 
Lithium Drifted diode, 1 mm thick, 
with a sensitive area of 1 cm 2 , and a 
weight of 0.232 g. Figure 1 contains 
a block diagram that shows the 
supporting electronics. Charge pulses 
generated from the energy deposited 
by incident particles are passed to the 
charge-sensitive amplifier, then to the 
discriminator and peak hold detector, 
and finally, they are converted from analog to 
digital signals by the 8-bit ADC (Analog-to- 
Digital Converter) and sorted into 256 channels, 
according to their voltage amplitudes. Two 
microprocessor units manage the data flow to 
and out of the memory devices. 


Table 1. Specifications and dynamic ranges 


the control of the program in the microproces- 
sors to collect dose, flux, and spectra until the 
memory is filled. At that time, the instrument is 
automatically turned “OFF” through the relay. 
The total time spent in the Working Mode is 480 
hours. The expected lifetime of the lithium 
batteries is about 1000 hours. 


Specifications 

Ranges 

Mass: (including batteries) 960 g 

Flux: 1 - 5000 part/cm 2 -s 

Size: 140 x 80 x 30 mm 3 

Energy Deposited: 0.16 - 13.8 

Temperature: -20-C to +45 Q C 

MeV 

Measurement Time: 30s to 30 min 

LET: 0.16-13.8 kev/p 


The specifications and dynamic ranges of the 
instrument are listed in Table 1 . The energy 
deposited or loss in silicon by the incident 
particles ranges from 0. 16 to 13.8 MeV. The 
measurement time per spectrum of the instru- 
ment can be varied from 30 s to 30 min. The 
LIULIN-3M operates in two modes, namely, 
“Working Mode” and “Data Transfer Mode.” In 
the Working Mode, the switch “ON” button 
connects the instrument via a jack to a polarized 
relay. In this mode the device is operating under 


In the Transfer Mode the device is switched 
“ON” after it is delivered to the principal investi- 
gator. This mode allows the transfer of the 
accumulated data in the flash memory to a 
personal computer (PC) under a parallel port 
communication link. During this mode the data 
are correlated to the real time of the flight. 

The major parameter of the instrument is the 
amplitude of the charge pulses generated in the 
solid-state detector by the incident particles. 






They are then passed to the charge-sensitive 
amplifier that converts them to corresponding 
voltage pulses. The sensitivity or conversion 
factor of this amplifier is 240 mV/MeV. The 
threshold setting of the discriminator is 38.7 
mV which establishes the lower data limit 
and consequently, the number of channels 
assigned to noise signals. Since there are 256 
channels in the Analog-to-Digital converter 
and the upper voltage limit is 3.3 V, the 
sensitivity per channel is 3.3/256 = 12.89 
mV/Ch and thus the first three channels are 
assigned to storing noise pulses. The depos- 
ited energy (Dep E) is given by equation 

Dep E(MeV) = £. (counts. * channel) * M ( 1 ) 
i = 4 - 256 

where the conversion factor M is given by: 

M={ 12.89^Vf240— 1=0.05371 ( 2 ) 

( Ch / { MeV J Ch 

and the corresponding dose rate per hour by: 

Dose Rate D R (pGy si /h) = K * Dep E (MeV)/t c (3) 

where t c is the collection time (in seconds) 
and the conversion factor K is: 

K = 2.486 in units of pGy si / h * MeV (4) 

m. RESULTS 

All data presented herein were collected on 
commercial flights every five minutes but the 
conditions of the measurements were not 
systematically controlled, particularly in 
terms of location and altitude. Figure 2 is a 
plot of dose rate versus minutes after the 
instrument was turned on for a flight from 
Washington, D.C., to Los Angeles, Calif. 

Figure 3 covers the trip from Los Angeles to 
Auckland, New Zealand, and Figure 4 shows 
the last leg of the flight from Auckland to 
Christchurch, New Zealand. 


Figures 5, 6, and 7 represent plots of measure- 
ments for flights from Washington, D.C., to 
Denver, Colo., Denver to Houston, Texas, and 
finally, Houston to Denver, then back to Wash- 
ington the same day, respectively. Figure 8 is a 
summary plot of these four flights. Figures 9 and 
10 show data taken on flights from Washington, 
D.C., to Paris, France, and also on the return 
flight. 

Figure 1 1 contains a plot of the instrument’s 
measurements of dose rate versus time after turn- 
on for the Sofia-New York flight, indicating the 
background at the airport and the three data 
points for which the altitude levels are given. 
Figure 1 2 contains a plot of the data showing a 
curve fit to the four available data points versus 
altitude. Figure 13 depicts the data for the flights 
from Washington to Budapest. The measure- 
ments obtained on the return flight segments are 
plotted in Figure 14. 

In most of the commercial flights, altitude 
information was made available by the airlines 
after the flights; however, for the trip from Sofia, 
Bulgaria, to New York, New York, altitude data 
were provided to members of the LIULIN team 
on board the plane by the pilot on three occa- 
sions, while on a recent flight from Washington, 
D.C., to Frankfurt, Germany, to Budapest, 
Hungary, flight information (latitude, longitude, 
altitude, speed, temperature) was continuously 
presented on TV monitors in the cabin. It should 
be noted, however, that for some flights altitude 
data could not be obtained. 

Table 2 presents the total doses measured for 
travel between the indicated cities. 

IV. DISCUSSION 

For all flights reported in Figures 2 to 14, the 
peak dose rates are a function of two variables: 
altitude and magnetic latitude. Thus, in Figure 2 
for the trip from Washington to Los Angeles the 
peak rate reached almost 0.2 jiGy/5-min, 
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whereas in Figure 3 for the trip from Los Ange- 
les to Auckland, the dose rate remained around 
0.07 fiGy/5-min for the major part of the flight 
over the Pacific Ocean, and only climbed to a 
peak of about .12 pGy/5-min near the end of the 
trip (last 2.5 hours). The dose rate for the short 
flight from Auckland to Christchurch, Figure 4, 
peaked at about 0.09 pGy/5-min. 
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Llulln-3M Commercial Flight Data 

Los Angeles, CA - Auckland, NZ (1 2/8/97) 
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Figure 4. 


In this context, particularly interesting are the 
measurements from the Los Angeles to Auckland 
flight, shown in Figure 3. The low dose rates 
would imply a correspondingly low altitude, 
perhaps 30000 ft, but the airline confirmed that 
this specific flight took place at an altitude of 
39000 ft. 

These measurements confirm that, as expected, 
flights toward lower magnetic latitudes, that is 
toward the magnetic equator, will experience 
smaller doses. The reason is that only a few very 
high energy galactic cosmic rays (several GeV/ 
nucleon) can penetrate to those regions of the 
globe on account of the Earth’s magnetic field 
which deflects most arriving particles with lower 
energies (rigidity concept). The low dose rates 
shown in Figure 3 are probably due to the shield- 
ing effects of the Earth’s magnetic field which 
compensated for the normally higher cosmic ray 
progeny distribution at the elevated altitude of 
this flight. 

The peak rates during the two flights covered in 
Figures 5 and 6, Washington-Denver-Houston, 
were about 0. 15 and 0.13 (i.Gy/5-min respec- 
tively. During the return trip, however, the peak 
rates achieved a surprisingly high level of 0.19 
pGy/5-min for the Houston-Denver segment and 
0.28 jiGy/5-min for the Denver- Washington 
segment, as shown in Figure 7. The summary of 
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these four flights is shown in Figure 8. The 
altitude of the aircraft, as given to us at a later 
date by the airline, was highest on the Houston- 
Denver (35000 ft) and Denver- Washington 
flights (37000) compared to the Washington- 
Denver (35000 ft) and Denver-Houston (33000 
ft) flight. It remains to be determined why the 
short duration Houston-Denver flight at 35000 ft 
experienced a higher dose rate (by ~ 20%) than 
the Washington-Denver flight at the same alti- 
tude but of much longer duration. 
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Liulin-3M Commercial Flight Data 

Denver, CO - Houston, TX 
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Figure 6. 



Liulin-3M Commercial Flight Data 

Houston, TX - Denver, CO - Washington, DC 

Channels 1 2 Through 256 
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Figure 8. 


Figures 9 and 10 show data for a trip to Paris, 
France, and return to Washington, D.C., respec- 
tively. The longer flight time of about 6 hours 
shows that the rate was about constant for the 
flight to Paris, approximately 0. 18 pGy/5-min, 
but on the return to Washington a constant rise in 
rate is observed, with a peak of about 0.23 (iGy/ 
5-min at the end of the trip as shown in Figure 
10. Since we have no reliable altitude informa- 
tion for these flights, it is difficult to interpret the 
data correctly. 
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The flight from Sofia to New York provided for 
the first time an opportunity to measure the 
variation of the instrument’s response in real- 
time as a function of altitude through the coop- 
eration of the pilot. Figure 1 1 shows a plot of 
dose rate versus time after turn-on. During the 
trip to New York the plane flew at three different 
altitudes of 31000, 35000, and 37000 ft. It can be 
seen that the detector measured three distinct 
dose rates. Theses three rates are plotted in 
Figure 12 versus altitudes, including the back- 
ground rate on the ground. A curve fit was made 
to the four data points and extrapolated to 40000 
ft where the dose rate approaches 0.5 pGy/5-min. 


In Figures 13 and 14 for the trip to Frankfurt- 
Budapest and return, again a similar pattern was 
observed as that seen in Figures 9 and 10 for the 
Washington-Paris and Paris-Washington flights, 
namely the trips east to Paris (Fig. 9) and to 
Frankfurt (Fig. 1 3) indicate a relatively constant 
level of exposure rate during the transatlantic 
voyage, whereas both return trips west show a 
gradual increase in the dose rate as the flights 
progress. Assuming that the flight corridor is the 
same for both directions, it is conjectured that 
the variation is altitude dependent, that is, lower 
altitudes going east and gradually higher alti- 
tudes going west, maybe in consideration of the 
jet stream (?). This will be investigated and 
resolved in a follow-on paper. 
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Liulin-3M Commercial Flight Data 
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Liulin-3M Commercial Flight Data 
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Figure 14. 


All figures show low “background” measure- 
ments on the ground before take-off and after 
landing, ranging from 0.017 to 0.033 (lGy/5-min. 
Upon close examination it is apparent that the 
data from low-altitude locations like Washington, 
Houston, and Los Angeles, with elevations less 
that 100 ft, have almost the same values as data 
from elevated locations as, for example, 
Snowmass, Colo., at over 9000 ft altitude (see 
Fig. 8), that is, average of about 0.0208 (iGy/5- 
min and 0.0225 pGy/5-min, respectively. We had 
expected to observe a larger difference at these 
two altitudes. The normal global average rate of 
radiation background levels, excluding high 
concentrations of naturally occurring radionu- 
clides such as radium, thorium, uranium, etc., 
has been reported to be about 0.008 to 0.017 


pGy/5-min at about 1 meter above ground and 
including cosmic ray contributions coming from 
above [5]. It has been established that the detec- 
tor orientation has no discernible impact on the 
measurements, in flight or on the ground. 

The total dose accumulation for each traveled 
flight segment has been calculated and is pre- 
sented in Table 2. While scheduled airlines are 
supposed to follow established flight corridors 
and altitudes, many times changes are made 
along the way because of safety, traffic, or 
weather conditions. Therefore, nominal flight 
plans, issued before the actual departure, which 
often may be delayed, are of no help in correctly 
analyzing and interpreting the measurements in 
terms of magnetic rigidity (cutoff latitude) and 
altitude. For this purpose, a “corrected” flight 
plan is required which indicates for every course 
change the exact location (geographic coordi- 
nates), altitude, zulu time, with several additional 
interim positional “fixes.” This type of data was 
not available for the initial flights discussed in 
the present paper but is now systematically being 
obtained for long-range global mapping purposes. 

As additional data are being accumulated over 
global distances with large longitudinal and 
latitudinal variations, a comparison to some 
predictions of cosmic ray intensities and their 
progeny may be possible, not only in terms of 
altitude but also as a function of rigidity and 
geomagnetic shielding affects. 

Therefore, in Table 2 we present only the total 
accumulated dose for each flight, the flight 
duration in minutes, and the average dose per 
minute. The highest exposure occurred on the 
Denver, Colo., to Washington, D.C., trip (3.84 x 
10 2 pG/m) and the lowest on the Auckland, New 
Zealand, to Christchurch, New Zealand, segment 
(1.1 Ox 1 O' 2 |iG/m), a factor of about 3.5 different. 
Other high doses were experienced on the 
Washington-Paris/Paris-Washington and the 
Washington-Frankfurt/Frankfurt-Washington 
flights. This is probably due to the fact that all of 
these flights traveled some of their time at high 
magnetic latitudes. 
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Table 2. Total Doses and Flight Duration 


Departure 

Airport 

Arrival 

Airport 

Figure 

Total Flight 
Dose (uGy) 

Flight 

Duration 

(minutes) 

Average Dose 
(pGy/m x 10 2 ) 

Washington, DC 

Los Angeles, CA 

2 

7.073 

260 

2.72 

Los Angeles, CA 

Auckland, NZ 

3 

10.483 

685 

1.53 

Auckland, NZ 

Christchurch. NZ 

4 

0.604 

55 

1.10 

Washington, DC 

Denver, CO 

5 

3.967 

185 

2.14 

Denver, CO 

Houston, TX 

6 

1.773 

90 

1.97 

Houston, TX 

Denver, CO 

7 

2.512 

95 

2.64 

Denver, CO 

Washington, DC 

7 

5.187 

135 

3.84 

Washington, DC 

Paris, France 

9 

12.018 

380 

3.16 

Paris, France 

Washington, DC 

10 

14.835 

435 

3.41 

Sofia, Bulgaria 

New York, NY 

11 

16.387 

555 

2.95 

Washington. DC 

Frankfurt. Germanv 

13 

12.739 

390 

3.27 


Budapest, Hungary 

13 

1.236 

60 

2.06 

Budapest, Hungary 

Frankfurt, Germany 

14 

1.219 

60 

2.03 

Frankfurt, Germany 

Washington, DC 

14 

14.721 

460 

3.20 


V. CONCLUSIONS 

The LIULIN-3M instrument is an effective 
radiation monitor for aviation or space applica- 
tion. It is characterized by small size, low power, 
light weight, and low cost. The instrument is 
currently being used to perform global surveys 
of radiation environment levels within aircraft 
and on high-altitude balloon missions over the 
South Pole (U.S. Antarctic Program). Measure- 
ments from a multitude of commercial airline 
flights conducted so far indicate a strong altitude 
and magnetic latitude dependence, as predicted 
by theory. The maximum dose rate encountered 
until now is about 0.3 |iGy per 5-minute interval. 

VI. REFERENCES 

[1] V.M. Petrov, V.S. Mahkmtov, N.A. Panova, 
V.A. Shurshakov, Ts. P. Dachev, J.V. Semkova, 
and Yu. P. Matviijchuk, “Peculiarities of the 
Solar Proton Events of 19 October 1989 and 
23 March 1991 According to the Measure- 
ments Onboard the Mir Space Station,” Adv. 
Space Res., Vol. 14, No. 10 pp (10) 645-(10) 
650, 1994. 

[2] Dachev, Ts., J. Semkova, V. Petrov, V. 

Redko, V. Bengin, T. Kostereva, J. Miller, L. 


Heilbronn, C. Zeitlin, “Analysis of the Pre- 
flight and Post-flight Calibration Procedures 
Performed on the LIULIN Space Radiation 
Dosimeter,” Acta Astronautica, 42, 375, 1998. 

[3] Dachev, Ts. P., J.V. Semkova, Ju. N. 
Matviichuck, B.T. Tomov, R.T. Koleva, P.T. 
Baynov, V.M. Petrov, V.M. Shurshakov, Yu. 
Ivanov, “Inner Magnetosphere Variations After 
Solar Proton Events. Observations on Mir 
Space Station in 1989-1994 Time Period,” 
Adv. Space Res., 22, No 4, 521, 1998. 

[4] Ts. P. Dachev, B.T. Tomov, G. Georgiev, P.G. 
Dimitrov, Yu. N. Matviichuk, V.M. Petrov, 

V.V. Shurshakov, E.G. Stassinopoulos, and J. 
Bart, “Analysis of the Calibration Results of 
the Engineering Model of LIULIN-3 Radiom- 
eter-Dosimeter,” Proceedings of the 3rd 
National Conference on Solar-Terrestrial 
Influences, 11-73, Sofia, Bulgaria, June 1998. 
(In Bulgaria) 

[5] P.F. Gustafson and S.S. Brar, “Gamma- 
Emitting Radionuclides in Soil,” The Natural 
Radiation Environment, Editors John A.S. 
Adams and Wayne H. Lowder, University of 
Chicago Press, 1964. 


7 









































































REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, Including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202 4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704*0188), Washington, DC 20503 


1. AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 

February 2002 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


4. TITLE AND SUBTITLE 

The LIULIN-3M Radiometer for Measuring Particle Doses in Space 
and on Aircraft 


6. AUTHOR(S) 

E.G. Stassinopoulos, C.A. Stauffer, T.P. Dachev, G J. Brucker, 
B.T. Tomov, and P.G. Dimitrov 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS (ES) 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 



8. PEFORMING ORGANIZATION 
REPORT NUMBER 

2002-01124-0 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS (ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 

TM — 2002-2 1 0003 


11. SUPPLEMENTARY NOTES 

C.A. Stauffer, Stinger Gaffarian Technologies, Greenbelt, Maryland; 

T.P. Dachev, B.T. Tomov, and P.G. Dimitrov, Bulgarian Academy of Sciences, Solar-Terrestrial Influences Lab, 
1113 Sofia, Bulgaria; G.J. Brucker, Radiation Effects Consultants, Inc., West Long Branch, NJ 07764. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified-Unlimited 
Subject Category: 35 

Report available from the NASA Center for AeroSpace Information, 

7121 Standard Drive, Hanover, MD 21076-1320. (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

This paper reports on the development of a compact radiation monitor/dosimeter, the LIULIN-3M, and 
on extended measurements conducted on the ground and on commercial aircraft on domestic and 
international flights. 


14. SUBJECT TERMS 


LIULIN-3M, radiation measurements, doses, radiation exposures. 


15. NUMBER OF PAGES 

7 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified UL 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39 18 
298-102 



















